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High-speed Terahertz communication systems has recently employed 
orthogonal frequency division multiplexing approach as it provides high 
spectral efficiency and avoids inter-symbol interference caused by dispersive 
channels. Such high-speed systems require extremely high-sampling 
time-interleaved analog-to-digital converters at the receiver. However, 
timing mismatch of  time-interleaved analog-to-digital converters 
significantly causes system performance degradation. In this paper, to avoid 
such performance degradation induced by timing mismatch, we theoretically 
determine maximum tolerable mismatch levels for orthogonal frequency 
division multiplexing communication systems. To obtain these levels, we 
first propose an analytical method to derive the bit error rate formula for 
quadrature and pulse amplitude modulations in Rayleigh fading channels, 
assuming binary reflected gray code (BRGC) mapping. Further, from the 
derived bit error rate (BER) expressions, we reveal a threshold of timing 
mismatch level for which error floors produced by the mismatch will be 
smaller than a given BER. Simulation results demonstrate that if we preserve 
mismatch level smaller than 25% of this obtained threshold, the BER 
performance degradation is smaller than 0.5 dB as compared to the case 
without timing mismatch. 
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1. INTRODUCTION 

In recent years, orthogonal frequency division multiplexing (OFDM) is a capable multi-carrier 
modulation technique that has been employed in terahertz (THz) communication systems because it saves 
transmission bandwidth and provides an ability to transmit high-speed data through dispersive 
channels [1]—[4]. However, at the receiver, regular analog-to-digital converters (ADCs) cannot operate at 
such extremely high sampling rates due to the limitation of the current semiconductor technology [5]. To 
overcome hardware constrains, ADCs with time-interleaved (TI) architecture are usually considered as a 
practical solution. In general, a time-interleaved (TI-ADC) consists of L parallel sub-ADCs, which 
simultaneously operate for the same period T,. In this fashion, the overall sampling rate of a TI-ADC is L 


i f 1 : er 
times higher than the sampling rate 7A of each sub-ADC. However, due to asymmetric manufacturing in a 
Ss 
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TI-ADC unit, the outcome data can be dramatically degraded by mismatches amongst sub-ADCs. One of the 
major problems of a TI-ADC is timing mismatch, which is induced by different sampling jitters of sub- 
ADCs. Over the past decades, timing mismatch’s impact and its calibrations were fully researched for single- 
carrier communication systems [6]—[17]. In recent past, when OFDM was employed in THz communication 
systems, this issue has been examined for multi-carrier communication systems, and it was shown that the 
timing mismatch affects OFDM system performance in a different behavior as compared to single-carrier 
communication systems [18]—[23]. To the best of our knowledge, the impact of timing mismatch on bit error 
rate (BER) performance in OFDM systems was investigated by simulation only in most studies, which can be 
very time consuming. This motivated us to analytically study OFDM-BER performance in the presence of 
TI-ADC timing mismatch. Further, in practice, the circuits-and-systems designing engineer is interested in 
determining tolerable timing mismatch levels which induce a negligible system performance degradation. We 
note that, until now, the identification of tolerable TI-ADC timing mismatch levels has never been 
investigated for OFDM systems. Hence, in this paper, we will determine tolerable timing mismatch levels of 
TI-ADCs for OFDM systems in Rayleigh channels. To this end, we derive approximate BER expressions 
based on the presumption that the interference induced by timing mismatch is Gaussian distributed. Most 
importantly, the derived BER expressions allow us to reveal a rule-of-thumb for identifying the maximum 
timing mismatch level, which can be condoned to assure that the performance degradation is smaller than 
0.5 dB as compared to the case without timing mismatch. 

The paper is constructed as follows. The BER expressions for square quadrature amplitude 
modulation (QAM) and pulse amplitude modulation (PAM) are obtained in section 3, assuming binary 
reflected gray code (BRGC) bit mapping [24]. In section 4, the accuracy of the obtained BER expressions is 
confirmed by comparing numerical results to Monte Carlo (MC) simulations. In section 5, we reveal a 
rule-of-thumb for tolerable timing mismatch levels producing an acceptable degradation for fixed timing 
mismatch and random timing mismatch. Section 6 will conclude the paper. 


2. OFDM SYSTEM MODEL WITH TI-ADC TIMING MISMATCH 

A block diagram of the considered OFDM system employing a TI-ADC at the receiver is shown in 
Figure 1. In the TI-ADC, the different sub-ADCs will have different timing error values, i.e., dt, denotes the 
timing error of the I“ sub-ADC. Let the vector X denote the input of an inverse discrete Fourier transform 
(IDFT) of size N. Hence, the vector X consists of N complexes-values symbols taken from M2-QAM 
constellation, i.e., X=(X0, X1, ..., XN-1)T, with each constellation symbol corresponding to a sequence of 
m = log2M bits, assuming the BRGC mapping rule [24]. Next, the cyclic prefix (CP) of length NCP is 
added to the IDFT output to protect the received data against inter-symbol interference (ISI) induced by 
dispersive channels. The output of the CP insertion unit sx is expressed as (1), 


aa. CK: 
Sk = eee ee = —Ncp, SNe 1,...,N— 1 (1) 
Channel 


Mappin OFDM Modulation Transmit 
pp'ng + CP insertion filter 
hy, 
AWGN noise 
Data OFDM Demodulation 
Detector + CP Removal 


Wr 
Figure 1. An OFDM system block diagram employs a TI-ADC 


Receive 
Filter 


Further, a digital-to-analog converter (DAC) at the transmitter will convert the discrete-time 
signal (1) to a continuous-time signal, which is then shaped by a transmit filter before passing through the 
channel hx with an additive white Gaussian noise (AWGN) noise wk. At the receiver, the received waveform 
is passed through a matched receive filter. Then, the 1“ sub-ADC samples the resulting signal at sampling 
time kLTs + ITs + dt,Ts, where 1=0, 1, .... L—1. In our analysis, we assume the TI-ADC to take a 
sufficiently high resolution in order that the quantization noise is ignored [19]. Moreover, in the most 
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practical cases, the timing errors dt, of each sub-ADCs change very slowly in time [19], and hence we can 
consider their values as constants during an OFDM symbol period. By employing the TI-ADC model 
affected by timing mismatch in [19], the TI-ADC output is written by (2), 


Te = SEs Enie Mm Sk—m(dt)) + Wo k = —Nep, Ncp + 1,..., N — 1, (2) 


where, Es denotes the symbol energy, L = mod(x, y) denotes the remainder after division of x by y, Amis the 
m™ tap of the sampled impulse response of the channel, w; are independent and identically distributed (i.i.d.) 
AWGN samples, i.e., wx ~N (0, No), and sx-m(dt)) is given by (3) 


Sk-m(dt) = E 1 y e?r m-at (3) 
Substituting (3) into (2), after appropriate re-arrangements, we obtain (4), 


ck 
Tk = Fs yn “3 X He -jan T fay + Wr (4) 


where, H, = ESP hy e? jane denote the frequency channel coefficients. Furthermore, the receiver will 
remove the CP and use a discreet Fourier transform (DFT) unit to convert the remaining N samples 


(from CP removal) to the frequency domain. The DFT output is expressed as (5). 


kn 
1 = -jizni 
Ry = g Eko ke Nn = 0,1, N — 1 (5) 
To further shorten (5), we consider the concept of a window function mg, which is defined as 7x is 


equal to 1 for k=0, 1, ..., N—1, and is equal to O for other cases. Its discrete-time Fourier transform (DTFT) 
II(f) is defined as (6): 


ICF) = NYF_.. sinc (N(f — k))e INF ~ Nsinc((Nf) ye IPN (6) 


where, (z)y = mod(z+ N/2, N) —N/2. Further, we substitute (4) into (5), and use the concept of the 
window function mg to change these summations to kE[—00, +00], and gL+l, gE[-o, +], /E[0,1, ..., L-1]. 
Therefore, we have (7): 


dt = 
L jon o)l qez o)qL 


a = EINA XH Li e N eI R iL onne N +W, (7) 


where, W, = 1/ VN X4 wpe” pay denotes AWGN samples in the frequency domain with W,,~N(0, No). 
Further, we consider the last summation of the first term in (7) as the DTFT of a time-shifted version of zz. 


. . . n-cjE é 
Hence, when we assess this version with f = To we obtain (8), 
fe (n-c)qL g et il 
—j20. j 2T — laL- (n=c)  i\ -jan= 
e” TOSS Taste” N= > Lino n (2 -: e74" (8) 


when we substitute (6) and (8) into (7), with proper arrangements, we have (9), 
n ~ EsTa XnHn + Yn = 0,1,...,N—1 (9) 
where (9) holds for large N, and Tn is given by (10), 


Tn = Did DTing(n — Pin) (10) 


zat} 
with DT, =-yicpe Pw to PTN, g(z) = ei", and Pin = mod (n- i7,N). Further, in (9), 


Yin is an ATE EE E term, which EEN as (11): 
Y= VY, +W, (11) 
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where, Y4 n is an inter-carrier interference (ICI) induced by TI-ADC timing mismatch, i.e., 


Yin = VE; Ving EN been DiveX cH gle = Pin) (12) 


à : 4 N = L-1 x 
Note that with integer ratios 7 we have Ta = DTo n and Yin = Xio PT; Pin Pin" pin: Moreover, assuming a 


perfect channel estimation, (9) becomes, 


(13) 


Finally, (13) are used to execute data bit sequence detection. 


3. BER DERIVATION BY USING GAUSSIAN APPROXIMATION 

Assuming that the symbols X, can be considered as statistically independent random variables with 
a uniformly distribution over a normalized-energy constellation, we have: i) Yin and W, are independent 
from each other and have zero mean and ii) ;,, is the summation of i.i.d. random variables. In the case of 
Rayleigh channels (i.e., H,~N(0,207)), each variable in the sum (12) is Gaussian distributed. As a result, 
Yin (12) has a Gaussian distribution. Hence, the variance oH, „Of Wn can be written as (14): 


Opn = = E; Dia tit 1 ce deyn D Ti, (DT, in, Ai ohg(c oo Di,n)g(c = Din) , (14) 


where, (z)* implies the conjugate of the complex number z. Therefore, the term Y, (11) is approximated as a 
Gaussian random variable, which has zero mean and variance per dimension. 


of = og, p, +2 (15) 


With the above approximation, for a given symbol Xn, timing mismatch dt, and a given channel H, 
we can obtain the conditional BER expression in an AWGN channel [21]. Further, in the most practical 
situations, not all N sub-carriers are data-modulated, i.e., some sub-carriers nearby the edges are not used to 
attain an efficient transition guard band [25]. Therefore, assuming only Nd sub-carriers for data transmission 
(Nd<N), we can obtain the BER for a given channel H, i.e., 


BER\y = —— Enuvy Cpa yerte ( Ynuv,y (Hn )) (16) 


NqamM 
where, the summation runs over n€{0,1, ..., Na-1}C{0, 1, ..., N-1}, uE{1, 2, ..., logoM}, vE{0, 1, ..., M-1} 
and ye {0, 1, ..., Guyv-1}. Moreover, in (16), Gu,» implies the number of decision boundaries left of the specific 
constellation symbol, which is given by (17): 


Guv = |(2v + 1) -2704+22 + 27t] (17) 


where, | z i indicates the largest integer smaller than z. Furthermore, erfc(z) in (16) is the complementary 


error function, and the argument Y; yy )(Hn) = (Spd — Bu yad )|Ta |, E >, with d the half minimum 


Euclidean distance [26], S, =2v+1-M and Byy=(2y-1) 2""*+!—-M. Finally, the pre-factors 
Quvy = (=1)l2"7?7™(Sv-Buy-=1)] in (16) take the values +1 or —1. It can be easily seen that because the BER 


for the J component or Q component of an M?-QAM constellation is similar, (16) also holds for an M-PAM 


modulation. From (16), we need to average the derived BER expression over the statistics of {|H|} in order 
_lHnl? 


to attain the BER for Rayleigh channels, i.e., BERpayteign = fs. BER X Hal e 2ofy d|H,,|. By using the 


result in [27], this integral can be further simplified as (18): 


BER = — 0, 


2 

OHE. 
n, u v,y |2 T = ) (18) 
rows OntoHEsOnuvy 
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where, On uvy =(S,d-B,,yd)|T,|. It can be obviously seen that (18) grants us to efficiently assess the system 
performance in the presence of TI-ADC timing mismatch. As a result, this will avoid a demand of very 
time-consuming MC simulations to evaluate system performance. 


4. BER EXPRESSION VALIDATION 

In this section, we will confirm the correctness of the obtained BER expressions in a Rayleigh 
channel by comparing (18) with MC simulations for a variety of constellations and mismatch levels. Further, 
in our simulations, we assume that all sub-carriers are data-modulated (i.e., N=Nd), and the sampling time Ts 
and the symbol energy E, are normalized. Moreover, L’ independent timing errors gre are uniformly 
generated over the interval [—1.1] [19] (see Table 1 with L = 8), which can be represented as the 100% 
timing mismatch level of a specific TI-ADC. Futhermore, when L<L’, only the first L timing errors of the L’ 
values in Table 1 are considered in simulations. Furthermore, we define the mismatch level by scaling 


dt)” in Table 1, i.e., for a mismatch level 2%, we have: dt7” = ai Finally, we consider the 


model of the channel impulse response in our simulations as (19): 


1 
hy, = Ke =C, k = 0,1,.....7-1 (19) 


where, T denotes the number of channel taps, C,;~N(0,1), and K is the normalization coefficient such that 
Vico Ef{|h,|7} = 1. In this case, it can be obviously seen that that Hp ~N (0,1). 

It can be seen from Figure 2 that a low mismatch level can even induce error floors in the BER 
performance. Further, Figure 2 shows that these error floors dramatically increase when either the mismatch 
level or the modulation order increases. Finally, Figures 2(a) and 2(b) reveals that the theoretical and 
simulated BER curves are in excellent agreement, which confirms the correctness of the obtained BER 
expressions. We also examined other parameter settings (results are not shown here) and identified the same 
excellent agreement between theory and simulation. This can be expected as we do not involve any 
approximation to derive the BER expressions for Rayleigh channels in our analysis. 


Table 1. Reference timing errors 


Parameter Reference values 
di}? [-0.92, 0.58, 0.49, 0.79, —0.17, 0.7, —0.44, 0.76] 


100% 
10" 
107 i 
492 
w 10 
a 
10* 
40-5 | —— No mismatch —#— No mismatch 
——  Mismatch-Analy. Expr. Mismatch-Analy. Expr. 
* — Mismatch-Simulation * — Mismatch-Simulation 
10% T T T T T T 
0 10 20 30 40 50 60 70 0 10 20 30 40 50 60 70 
E/N, (dB) E,/Ny (dB) 
(a) (b) 


Figure 2. BER system performance for a Rayleigh channel: (a) L=8, N=2048 and 1% mismatch for different 
modulation types and (b) 8-PAM, N=2048 and L=7 with different mismatch levels 


5. AN IDENTIFICATION OF TOLERABLE TIMING MISMATCH LEVELS 

In this section, the tolerable mismatch level inducing a negligible performance degradation will be 
identified. To this end, we first derive a condition for the maximum timing mismatch level for which error 
floors induced by timing mismatch will be smaller a given BER value, i.e., BER,. Next, from this derived 
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condition, we will reveal a rule-of-thumb to determine tolerable timing mismatch levels at which the system 
performance degradation is not greater than 0.5 dB at a given BER value as compared to the mismatch-free 
case. Section 4 illustrates that the timing mismatch consistently induces the error floor at high SNRs. Further, 
we interpret the maximum timing mismatch level yg% for which error floors are lower than BER;. Moreover, 
let us define the tolerable timing mismatch level fR% as the maximum level for which, at BER;, the BER 
performance degradation is not greater than 0.5 dB as compared to the case without mismatch. Typically, we 
have 7rp% < yr%. Considering (18), to assure error floors to be lower than BER, at SNRs, we consider the 
dominating BER terms only, i.e., 


_lHnl? 
ZENA max fy” erfe(¥(y%)) x “tle 2h dlH,| < BER, (20) 
N uwy `O on 
M-1 % [__&s 
where, B = MGA? Y (y%) = (S,d = Big) |T EA 208, „0% and 
03, , (2%) = 
E; Èi ize{o,1 a L-1}\0 È ceto, sa N-1},c#n DTZ% (DTZ%) onf (c z Pin Fle = Pin) (21) 


by using the results in [27], (20) can be simplified as (22), 


2 
(Syd-By,y a) |T %| oR Es 


5 N 
Ed max 41 -— ae te 5 BER: (22) 
Uv, aĝ, n(V%) +O} Es (Spd-By yd) ee 
or equivalently 
g a 
=e 0b, ,(V%) 5 
Ng mins | < a + of E,) = Jr (23) 
uwy | \(Syd-Bu yd) [ner 


where Jp = N(ofE,)~*/? (1 — e, From (23), it shows that we minimize the left hand side over (u, v, y) 


which means minimizing S,d—B,,,d corresponding to (u, v, y). Because S,d—B,,)d denotes a distance between 
constellation symbols and decision boundaries, we have min{S,d—B,,yd}.,1,y=d. Therefore, (23) becomes (24), 


2 


2 = 
= Ob, (YR) 
neo [se + cit. 2 Jr (24) 


2|p ÝR% 
d ees 


Taking into account (24), the threshold yp%, at which error floors produced by the mismatch is lower than 
BER,, will be next explored, i.e., if we select are in Table 1 as y% x dtt% considering y% < Yr% and 
BERz,/n +0 S BER;. To find the threshold yg%, we examine (24) with an equality. From this obtained 
threshold level yp%, the maximum mismatch level 7%, which induces a performance degradation not 
greater than 0.5 dB as compared to the case without mismatch, will be determined. Figure 3(a) demonstrates 
that in all considered cases, if FR% < 0.25yp%, there is an acceptable performance degradation (<0.5 dB) at 
BER, = 107°. Further, Figure 3(a) also illustrates that if the mismatch level increases to yp%, there is a 
noticeable performance degradation. We also investigated for other number of sub-ADCs, and other 
modulation types and orders, and found the similar results. As a result, we can conclude a rule-of-thumb for 
tolerable timing mismatch levels as (25). 


Fr% < 0.257 p% (25) 


Until now, we obtained the results with fix values dt, only. Nevertheless, in the most practical 
situations, the timing error values dt, are random. Therefore, we will assess the average BER performance 
with random values dt,. This can be attained by averaging the BER system performance over the larger 
number of different TI-ADCs. In the following results, the timing error values dt, are uniformly chosen over 
[-z/100, z/100], 0<z<100, where z determines the mismatch level in percentage. Further, in each simulation, 
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we will choose different sets of L timing error values dtl. Moreover, the results in Figure 3(b) are obtained 
by employing 1010 TI-ADC realizations, BER, = 10~°, and different mismatch levels. It can be seen from 
Figure 3(b) that the theoretical and simulated BER curves do match well in a Rayleigh channel. Most 
importantly, Figure 3(b) demonstrates that at BER;, if the mismatch level is equal to the tolerable threshold, 
i.e., O.25yp%, where yp% is found by evaluating (24) with an equality, the performance degradation is 
invisible as compared to the mismatch-free case. On the other hand, there is a noticeable performance 
degradation if we select the mismatch level larger than 0.5yg%. Hence, we conclude that the obtained 
rule-of-thumb of tolerable mismatch levels is valid for Rayleigh fading channels. 


—— No mismatch 
© 0257% 
A- 05y% 

-A 7% 


— No mismatch 
© 0.25y,% - Analy. Expr 
% =0.25y,% - Simulation 

= - 0.5y,% - Analy. Expr 

A 0.5y,% - Simulation 

Ya% - Analy. Expr 

m y% - Simulation 


60 70 80 90 100 110 120 130 
E/N, (dB) E,/N, (dB) 


(a) (b) 


Figure 3. BER system performance for a Rayleigh channel with N=2048 and L=8: (a) fixed timing error 
values dt; and (b) 10!° sets of TI-ADC realizations 


6. CONCLUSION 

In this paper, the maximum threshold of TI-ADC timing mismatch level, for which error floors 
induced by the mismatch is not greater than a considered BER value, has been theoretically identified for 
OFDM systems in Rayleigh channels. Further, we have revealed that when the mismatch level is selected to 
be smaller than 25% of the reliably recommended threshold, the BER performance degradation is not greater 
than 0.5 dB as compared to the case without mismatch. This tolerable timing mismatch level can serve as a 
significant instruction for circuits-and-systems designing engineers to calibrate the TI-ADC timing mismatch 
through digital signal processing or hardware circuits. 
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